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Some high-energy experiments have suffered from synchrotron-radiation background. As a mea- 
sure, we have deveioped a new cakuiation method of synchrotron radiation based on a real beam 
orbit, aiming at quantitative estimations and the construction of a possible alarm system for the 
background. An explanation and a demonstration of our method are given. 



INTRODUCTION 



A. Fitting procedure 



At the initial stage of the KEK B-factory (KEKB) ex- 
periment synchrotron radiation (SR) from the cir- 
culating 8-GeV electron beam caused serious damage to 
the Belle silicon vertex detector (SVD) Q, which was 
located close to the interaction point (IP). Some of the 
preamplifiers mounted on the inner-most SVD layer died 
in about ten days just after start of the gain drop. Other 
high-energy experiments have also suffered from SR back- 
ground. 

We have developed a new method to calculate SR 
wattages based on a real beam orbit, aiming at quan- 
titative estimations and the construction of a possible 
alarm system for SR background. Here, the real orbit is 
obtained by fitting measurements of beam-position mon- 
itors (BPMs). In the following sections, our method is 
explained together with some illustrations using the 8- 
GeV electron beams of KEKB High Energy Ring (HER). 
Finally, we present a simulation of the orbit and the gain 
drop at the time of the SVD gain-drop accident as a 
bench mark. 



II. 



ORBIT CALCULATION 



Beam orbits are calculated by a linear approximation 
using transfer matrices, which are shown in Appendix lAl 
For numerical calculations, we coded a dedicated com- 
puter program in Fortran, where the correctness of the 
orbit calculation was checked by comparing with a cal- 
culation using SAD Furthermore, we confirmed that 
there was no significant effect from round-off errors by 
comparing the results of double- and quadruple-precision 
calculations. 
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Orbit fitting is performed based on the optics, field 
strengths of the magnets, and BPM measurements. We 
require minimization of x 2 > defined as 
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where xf PM) (yf PM >) 
beam position measured with BPM:j, X 



indicates the horizontal (vertical) 
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(orbit) \ 



the horizontal (vertical) position of the calculated orbit 
at BPM:j, and o~j is the BPM:j resolution. There are 
four floating parameters, which are the orbit positions at 
the entrance and exit of the relevant accelerator section: 
x (»n) ) y(»n) )J| r( ut) j y(o«t) < Minimization of x 2 is performed 

using the computer program MINUIT 0, 0] . An example of 
the orbit-fitting results for KEKB HER is shown in Fig. 
Here, we define the vertical and horizontal \ 2 to estimate the 
goodness of the fit for each direction separately: 
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which lead to the following relations: 
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indicates the number of degrees of 



freedom related to the horizontal (vertical) direction, and 



N df for both: N, 
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In the example 



shown in Fig. Q the goodness of the fit is Xx /Ndf(x) = 10.5 
and Xy/NdfCY) = 40.3 for a conservative BPM resolution of 
100 /im. As is often the case, \ 2 i s too large, and also the orbit 
at the IP is unrealistically deflected in the negative direction 
by about 3 mm. This is due to the fact that the absolute po- 
sitions of magnets and BPMs are not known sufficiently well, 
leading to a need for the offset corrections. 
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FIG. 1: Example of the orbit-fitting results for KEKB HER 
around the IP (s = AX = AY = 0) without any correc- 
tion. The left (right) figure shows the horizontal (vertical) 
fitted orbit together with the BPM measurements. The s 
axis is taken to be along the orbit. The horizontal line with 
AX = AY = corresponds to the design orbit. The HER 
beam flows from the left to the right. Xx/^dftx) = 10.5 and 
X 2 Y/N df(Y) = 40.3. 



B. Offset corrections 



We modify the \ 2 formula Q as follows: 
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where AXj^ 17 * - 1 (Ay^ UAD ') indicates the horizontal (verti- 
cal) offset, with respect to the design position, of the real field 
center of the quadrupole magnet (QUAD) on which BPM:j 



is mounted, and AXf r ' (AY, 



CPU , X ^(BPM) ) horizontal (verti _ 



.(QUAD) ^y-(QUAD) 



cal) offset of the real origin of the BPM:J measurement with 
respect to the real field center of QUAD:j. The index i runs 
for different orbits in the case of a global fitting, which is 
explained in the next paragraph. 

In order to obtain the offset sizes ( AXj . _ . 

AXj BPM ' and AY^ BPM '), we perform a global fitting for dif- 
ferent orbits within a time period in which magnets are not 
expected to have moved. In case of KEKB, we make vari- 
ous measurements in machine studies using six different or- 
bits in each of the horizontal and vertical directions, which 
are produced by making a small kick at one of the six steer- 
ing magnets. We fit those six orbits simultaneously with not 
only the usual floating parameters (X^"^ , Y^ m ' , Xf° , Y^°"'' 
(i — 1,2,3,4,5,6)), but also the offsets in formula ©. We 
furthermore introduce an additional floating parameter: the 
correction factor of the field strength of the bending magnet 
(BEND), a, which is defined as 8 corr — (1 + a) x 9 org , where 
6 org and 6 corr are the kick angles before and after the correc- 
tion, respectively. The BPM and QUAD offsets as well as the 
correction factors of the BENDs are common among the six 
orbits. It is clear that the simple fitting procedure leads to 
failure because there are too many floating parameters (more 
than one hundred in our example). 



FIG. 2: Example of the orbit-fitting results after the off- 
set corrections. The left (right) figure shows the horizontal 
(vertical) fitted orbit together with the BPM measurements. 



Xx/Ndf(x) 



3.6 and Xy/^/fy) = 0.6. 



Therefore, we first categorize the offset and correction- 
factor parameters to be floating into three sets: the BPM- 
related, QUAD-related, and BEND-related ones. We choose 
zero or one parameter to be floating in each category, and the 
other parameters are fixed at zero. Then, a global fitting is 
performed for the combination of the choice. We try all of the 
possible combinations (about 10 million in our example), and 
a large number of results are filtered, requiring the following 
criteria: 

• the minimization of x' 2 converges with a positive- 
definite Hessian matrix, 

. |AX] QUAD) | < 10mm and |AY/ QUAD) | < 10mm, 

• |AX, (BPM) | < 10 mm and \AY} BPM) \ < 10 mm, 

• |a| < 0.2, 

• the vertical and horizontal displacements of the fitted 
orbit at the IP with respect to the nominal IP are within 
1 mm. 

Then, the best (lowest) x 2 sample is selected as the final one. 
Figure |21 shows the fitted orbit after the offset corrections. 
With the corrections, the x 2 is small, and the orbit passes 
close to the BPM measurements and the nominal IP. Choosing 
two offsets to be floating in each category, we obtained almost 
the same results. 

This method is to find a solution with a minimum modifi- 
cation within the scope of the numerical approach. 



III. 



WATTAGE CALCULATION 



SR wattages are calculated, based on the fitted real orbit, 
using the following analytical formula |6j: 
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where a indicates the fine-structure constant, / the beam cur- 
rent, e the charge of the beam particle, 7 the Lorentz factor, 
ui the angular frequency, ui c the critical angular frequency, 4> 
the bending angle, ip the angle to describe the deviation from 
the bending plane, £ = 'yip, r\ = (1/2)(uj/lu c ) (l + £, 2 ) 3 ^ 2 , h 
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FIG. 3: Integration variables: <p and ip used in the wattage- \ t | | 

calculation formula J7J. The circle indicates a circulating -20 -10 10 20 

beam. X [mm] 



Planck constant (reduced), and K 2 /3 and .K1/3 are modified 
Bessel functions. Figure |H] explains the two variables: <p and 
ip. We regard A(j> as being an infinitesimal quantity, and the 
integral form for formula J7J is approximated by the following 
three-dimensional integration: 
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where the s axis is taken to be along the beam orbit, and 
p(s) is the radius of curvature at s. A wattage distribution in 
a plane can be obtained by taking the mapping of (s, ip) — > 
(X, Y), where (X, Y) is a point of intersection of the plane 
and the line with the angle ip with respect to the tangent of 
the orbit at s. Integration on s and ip corresponds to surface 
integration of the Poynting vector in the plane. We perform 
numerical integrations of formula © using the Monte-Carlo 
(MC) integration program BASES Q. 

The implementation of the wattage calculation in our pro- 
gram was checked by comparing with the numerical calcula- 
tion using SAD |8|. In the SAD calculation, the exact formula 
in electrodynamics is used. 

It is important to consider any beam-size effects in the SR 
simulation, especially for the SR from QUADs. The effects 
are included by adding four additional integration variables 
(X, X ', Y, Y'), which indicate a phase-space point of the orbit 
at the IP. MC integration is performed together with the fol- 
lowing Gaussian weight function to describe the beam profile: 



FIG. 4: Example of the wattage distributions at the IP in the 
plane perpendicular to the beam axis. The X axis is taken to 
be in the horizontal direction, pointing to the outside of the 
KEKB ring, and Y in the vertical direction. The SR comes 
from the vertical BEND located about 25 m upstream from 
the IP. The circle in the figure corresponds to the inner surface 
of the IP chamber. 



X'-Y' space. The beam sizes are calculated according to the 
following formulas: 



°x> = \/ex/y<3 
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°Y' = sj £y//3 y , 
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assuming a* = 0, where a* is one of the Courant-Snyder 
parameters at the IP in the notation: "/x 2 + 2axx' + fix' 2 = 
e. In the above formulas, ex (ey) indicates the horizontal 
(vertical) emittance, and X (/3 Y ) is the horizontal (vertical) 
beta function at the IP. The orbits in the relevant magnets 
are calculated using the inverse transfer matrices from the 
IP. This implementation in our program has been checked by 
comparing the results with no beam size and with a beam size 
in the asymptotic behavior: ex — > and ey — > 0. 

In Fig. 0] an example of the wattage distributions at the IP 
is shown in the plane perpendicular to the beam axis. The SR 
in this example was a dominant source in the SVD gain-drop 
accident. 
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where (Xfit,X'f it ,Yfit,Yf it ) indicates the phase-space point 
of the fitted orbit at the IP, ax (oy) is the horizontal (ver- 
tical) beam size, and a x > and oy> indicate the beam size in 



IV. 



BENCH-MARK TEST 



In order to demonstrate the validity of our method, we 
adopt the SVD gain-drop accident as a bench mark. For 
radiation-dose calculations, we made a detector simulation 
based on EGS4 IKl, which was improved for low-energy pho- 
tons [TH IT21 173. Il4l Il5| . In this detector simulation, the ex- 
act geometry and materials of the relevant SVD components 
have been installed according to engineering drawings. SR 
photons, generated with the general-purpose event generator 
SPRING [3, are fed into the detector simulation, and energy 
deposits in the SVD are calculated. Photons are simulated 
down to 1 keV, and electrons to 20-keV kinematic energy. In 
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FIG. 5: SVD gain as a function of time around the SVD 
gain-drop accident, together with the two periods used in the 
bench-mark test. The ladder numbers in this figure corre- 
spond to the position numbers in Tables Q] [H] and Fig. Q 
The gain-drop plots are extracted from [l7|l . 



order to transfer the radiation dose into the gain drop, we use 
the result on the irradiation test to investigate the radiation 
hardness of the SVD preamplifiers flrj . 

We make two different sets of calculations for a more care- 
ful demonstration. One of them is on the radiation dose be- 
fore starting the gain drop (period I), and the other is after 
the start (period II), as shown in Fig. [S] During each of the 
two periods, the beam orbit was almost unchanged. Table Q 
gives the results of a radiation-dose calculation for period 
I. Obtaining information on the time of the period (4641 
minutes), and on the average beam current (4.6 mA), the 
gain drop is estimated to be around or less than 1 %, 
which is consistent with the measurements (< 1%). Af- 
ter starting the gain drop, the estimated doses are much 
higher in period II than those in period I, as shown in Ta- 
ble[nl where the time of the period is 8071 minutes, and 
the average beam current is 6.7 mA. What changed sig- 
nificantly are the kick angle of the vertical BEND, which 
made a dominant contribution to the radiation dose in 
period II, and the orbit around the BEND, as shown 
in Fig. El The bending angle changed from 0.53 mrad 
with 1.7-keV critical energy (E cr i t ) to 0.66 mrad with 
E cr n = 2.2 kcV in the same direction. Figure shows 
a comparison between the estimated and measured gain- 
drop values. We have obtained good agreements not only 
on the absolute size, but also for the azimuthal distribu- 
tion. 



V. CONCLUSIONS AND FUTURE 
PROSPECTS 

We have developed a new method to calculate SR 
wattages, or radiation doses, based on a real beam or- 
bit with reasonable offset corrections. Here, the orbit is 
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FIG. 6: Fitted beam orbits around the SVD gain-drop acci- 
dent. The dotted line indicates the orbit in period I and the 
solid line in period II. The vertical arrow points to the BEND 
which made a dominant contribution to the SVD gain drop 
in period II. 
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FIG. 7: Comparison between the estimated and measured 
gain-drop values for period II. 



obtained by fitting BPM measurements. SR wattages are 
calculated based on the fitted orbit with the MC integra- 
tion of the analytical formula. We have coded a ded- 
icated computer program, where the correctness of the 
orbit and wattage calculations in our program has been 
confirmed by comparing with the results using SAD. Fi- 
nally, we have successfully reproduced the SVD gain-drop 
accident quantitatively, so that the practicability of our 
method has been established. 

This method can be applied to the construction of an 
online alarm system for SR background, leading to prior 
notice for detector damage and possibly more flexibility 
in accelerator operation. With this technique, a warning 
can be issued when the orbit is approaching dangerous 
areas, and furthermore the expected amount of SR back- 
ground can be estimated before making certain changes 
in the orbit. 

We have made a tentative version of the online program 
for KEKB operation. It takes about twenty seconds per 
cycle (getting the magnet information and BPM mea- 
surements, fitting the orbit, and MC integration of SR 
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Position number 


1 


2 3 4 


5 


6 7 





Dose [Rad/min/A] 
Integrated dose [kRad] 
Estimated gain drop [%] 


2.1±0.9 
0.06 
0.02 


1.2±0.7 23±3 41±4 
0.04 0.53 0.93 
0.01 0.17 0.30 


33±4 
0.75 
0.24 


33±4 25±3 
0.74 0.57 
0.23 0.18 


3.5±1.2 
0.10 
0.03 



TABLE I: Results of a radiation-dose calculation based on our method for the period before starting the gain drop (period 
I). The position numbers correspond to the azimuthal coordinate with respect to the beam axis, and the numbers 2, 4, 6, 
indicate the top, inside, bottom and outside of the KEKB ring, respectively. The errors from the MC statistics are shown 
beside the dose values. 



Position number 


1 


2 


3 


4 


5 


6 


7 





Dose [10^ Rad/min/A] 


1.5±0.4 


1.7±0.5 


10±1 


17±1 


24±2 


19±1 


19±1 


7.2±0.9 


Integrated dose [kRad] 


8.0 


9.2 


55 


92 


131 


102 


102 


39 


Estimated gain drop [% 


2.6 


2.9 


17.6 


29.4 


41.9 


32.6 


32.6 


12.5 



TABLE II: Results of a radiation-dose calculation for the period after starting the gain drop (period II). 



wattages from all the magnets in the HER straight sec- 
tion near the IP) using the workstation with a 500MHz 
CPU. If we use a recent faster computer, the turnaround 
time could be comparable to the time of BPM measure- 
ments (several seconds). 



QUAD (horizontal focusing) : 
fc -1 sin(/cL) 



/ cos(fcL) 
— k sin(kL) 
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cosh(fcL) 
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APPENDIX A: TRANSFER MATRICES 



The exact forms of the transfer matrices used in the 
orbit calculation are shown in this appendix. Here, L 
indicates the e ffective length of the component, k is de- 
fined as k = \J (B /b)/(Bp), where B is the magnetic 
field strength at radius b, and Bp is the magnetic rigidity 
of the central reference trajectory. 6 X (9 y ) is the horizon- 
tal (vertical) kick angle in the BEND. R rot (a) indicates 
the rotation matrix for the rotation angle, a. The ele- 
ments with a blank mean zero. 



Orbit : X(s) = 



' X'(s) * 
Y(s) 
Y'(s) 

\ 1 



(Al) 



Z 1 



Drift space : R-drift (^) 



(A2) 



f 1 



BEND : R bond = Rdrift(-) x 



X 



xR drlft (-) (A4) 



—J 



Skew QUAD : R BkowQ = R rot (-45°) x R quad x R rot (+45°) (A5) 
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